Cumulene Molecular Wire Conductance from First Principles 
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We present first principles calculations of current-voltage characteristics (IVC) and conductance 
of Au(lll):S2-cumulene-S2:Au(lll) molecular wire junctions with realistic contacts. The trans- 
port properties are calculated using full self-consistent ab initio nonequilibrium Green's function 
density-functional theory methods under external bias. The conductance of the cumulene wires 
shows oscillatory behavior depending on the number of carbon atoms (double bonds). Among all 
conjugated oligomers, we find that cumulene wires with odd number of carbon atoms yield the high- 
est conductance with metallic-like ballistic transport behavior. The reason is the high density of 
states in broad lowest unoccupied molecular orbital levels spanning the Fermi level of the electrodes. 
The transmission spectrum and the conductance depend only weakly on applied bias, and the IVC 
is nearly linear over a bias region of ±1V. Cumulene wires are therefore potential candidates for 
metallic connections in nanoelectronic applications. 

PACS numbers: 85.65.-(-h, 73.63.-b, 31.15.Ne, 03.65.Yz 



I. INTRODUCTION 

Molecular nanowires and atomic chains are actively ex- 
plored as novel conductors for nanoelectronic devices. 
Early on, cumulenes carbon chains were proposed as 
ideal molecular wiresi, and the interest in the struc- 
ture persist^>^. Subsequent work investigated the trans- 
port properties of cumulene wires directly connected to 
jellium^'^ or atomically structured^'^ electrodes via car- 
bon double bonds. For both types of electrodes, the con- 
ductance was found to vary in an oscillatory manner with 
the number of carbon atoms in the finite chain. Neverthe- 
less, the transmission not only depends on the properties 
of the wire, but also on the contacts to the electrodes. 
For example it was understood^"^ that the amount of 
charge transferred from the electrodes to the molecular 
wire chemically bonded to the metal electrodes, deter- 
mines the character of the conductance and its length de- 
pendence. The charge transfer from the electrodes to the 
cumulenes provides doping of wires^. Since reliable meth- 
ods of purification of cumulenes have been establishedSiS, 
we reconsider the transport properties of these molecules 
using a realistic metal-molecule interface model. 

In this paper we report full self-consistent ab ini- 
tio calculation of electron transport through the cumu- 
lene molecular wires =(C=)„ with n=4-9 connected to 
Au(lll) surfaces via thiolate bonds. Unlike the previous 
studies of the cumulenes^"—, where the cumulative double 
bond order was supported in the metalMcumulene form, 
with the metal atom at the electrode surface incorporated 
into the carbon double-bond framework, we here consider 
the full atomic structure of 'simple' organic cumulenes of 
type R2C=(C=)n-2CR2, forming metallOcumulene com- 
plex with the metal surface. Metallic-like conductance 
of another form of carbon chains with alternating sin- 



gle and triple bonds, polyynes, was recently studied by 
Crljen and Baranovic^*' using thiloate contacts; however, 
such a realistic model for the metallOcumulene wires was 
never investigated, because the molecule itself was con- 
sidered to be difficult to isolate in macroscopic quantity. 



II. COMPUTATIONAL METHOD 

The calculations are performed using density func- 
tional theory-based nonequilibrium Green's function 
(NEGF) transport theory as implemented in the Tran- 
SIESTA simulation package^. We used local-density ap- 
proximation (LDA) for the exchange correlation poten- 
tial. For carbynes, LDA performed better than gen- 
eralized gradient approximation (GGA) for calculating 
the structuresii. Using LDA for relaxation and recon- 
struction on Au(lll) surface has been extensively dis- 
cussed in previous literature^^. For the transport prop- 
erties one can expect the LDA to perform quite well due 
to the metallic character of the junctioniS,. Core elec- 
trons were modeled with Troullier-Martinsi^ soft norm- 
conserving pseudopotentials and the valence electrons 
were expanded in a basis set of local orbitalsi^. The 
transmission spectrum, which gives the probability for 
electrons with incident energy E to be transferred from 
the left to the right electrode, is calculated from^^ 

T{E, V) = tr [rR{E, V)Gc{E, V)rL{E, V)Gi{E, V)] (1) 

where Gc is the Green's function of the central region 
and Ti^/n is the coupling matrix of Left/Right electrode, 
respectively. The integration of the transmission spec- 
trum yields the electric current 
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T{E,V){f{E~tJiL)~ f{E-tin)}dE (2) 
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where fiL = —V/2 {fin = V/2) is the chemical potential 
of the left(right) electrode. 




of the cumulene systems with 4-7 carbon atoms in the bias 
FIG. 1: The zero bias electron transmission spectra of cumu- region from -2 to 2 V 
lene systems consisting of n=4-9 carbon atoms. 



In our calculations, free thiol-capped cumulenes were 
first optimized using the SIESTA packageiS. The bond 
lengths and geometrical structures obtained in all the 
performed calculations are similar to those of unsubsti- 
tuted cumulenesiZ.. The minimum-energy conformation 
of even-number wires, is planar, whereas the terminal 
groups of odd-number wires are mutually perpendicu- 
lar. A system of a cumulene wire coupled to Au(lll)- 
(3 X 3) surface electrodes via thiolate bonds was fully 
optimized^*. The adsorption geometry has been deter- 
mined by positioning thiolate capped cumulenes above 
four adsorption sites, on-top, bridge, fee and hep on Au 
surfaces at a favorable Au-S bonding distancei^. During 
relaxation sulfur atoms move to the bridge sites. Then 
the metal-molecule-metal sandwich was constructed, and 
fully optimized. The Au-S equilibrium distance is 2.20 A. 

To establish a firm ground concerning cumulene ad- 
sorption geometry, cumulene molecules on Au(lll) for 
all n considered were optimized by using the VASP— 
simulation package with LDA potential, following the 
computational procedure described in Ref. [20l : then the 
metal-molecule-metal sandwiches were constructed. We 
obtained identical geometries, as compared to SIESTA 
calculations, with marginal differences noticed only for 
the longest molecules. It is instructive to compare the 
sulfur adsorption geometry obtained in this work to the 
one obtained in Ref. [2^ The previously calculated ten- 
dency was for sulfur atoms to move from in-hollow to 
bridge and then to on-top position on Au(lll) when the 
sulfur density and/or hydrogenation density increases. 
This perfectly matches the present findings where two 
sulfur atoms brought together with bonds to a common 
carbon atom (but located relatively far from other Au:S 
assemblies) slide toward bridge positions from any initial 
adsorption geometry. For the sandwiched molecule, the 
bonding energy was found to be 1.7 eV per S atom for 
even n and 1.86 eV per S atom for odd n. For ri=8,9 
these energies becomes 0.1 eV smaller. 



III. RESULTS AND DISCUSSION 

Fig.[T]shows the zero bias electron transmission spectra 
(T) for the studied metal-molecule-metal systems, while 
the corresponding IVCs and difii'erential conductances are 
shown in Fig. [2] in the bias region from -2 to 2 V. The 
zero bias conductance in Fig. [3] clearly demonstrates the 
oscillating behavior of the conductance of odd/even-n cu- 
mulene chains. 

To examine the results in detail, we concentrate on the 
tetraene and pentaene chains. Fig. |4] shows transmission 
spectra and molecule projected self-consistent Hamilto- 
nian (MPSH) eigenvalues of the cumulene systems (a) 
tetraene (b) pentaene at zero bias. Atoms in the central 
region included in the self-consistent cycle are shown in 
the lower inset. The MPSH eigenvalues are in fair agree- 
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FIG. 3: Zero-bias transmission T{E — 0, n) through the car- 
bon chains. 
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FIG. 4: Transmission spectra and molecule projected self-consistent Hamiltonian (MPSH) of cumulene systems (a) tetraene 
(b) pentane at zero bias, as compared to energy levels of corresponding radicals. Atoms in central region included in the 
self-consistent cycle are shown in the inset. The broad resonances are marked as peak A and B at the first peak below and 
above the Fermi level, respectively. 



ment with the energies at the transmission resonance 
peaks. The main features of transmission spectra are 
the broad resonances below and above the Fermi level, 
marked as "peak A" and "peak B", respectively. Com- 
paring the MPSH states with molecular orbitals (we use 
deprotected radical for clarity) we conclude that the den- 
sity of states contributing to the peaks A and B arises 
from the unoccupied levels, and that the LUMO level 
dominates transmission channels. Undoubtedly, the high 
value of the conductance of cumulenes is to be expected 
because of the delocalized nature of the LUMO. 

Oscillatory behavior of conductance with the number 
of carbon atoms in the chain is evident from Fig. |31 with 
odd-n carbon atom wires giving a higher conductance 
than even-n ones. All cumulenes demonstrate high con- 
ductance, higher than other conjugated polymers, and 
for the odd-n cumulene wires the conductance is even 
higher than polyyues^^ (see Fig. [S]). The high value of 
the conductance results from the high density of states at 
the Fermi level, with the position of the LUMO aligned 
with Fermi level, as clearly seen in Fig. |4l While the 
conductance of polyynes slowly decreases with increas- 
ing wire length, the conductance of odd-n cumulenes does 
not show any pronounced dependence on the molecular 
length, which is a clear signature of ballistic conduc- 
tance. Because of ballistic transport, backscattering in 
molecular wires should be negligible, with little energy 
dissipation, with finite resistance attributed to the effect 
of backscattering from the contact region (constriction 
and interface structure). Previous work^ on cumulene 
systems with jellium electrodes has found Schottky-like 
barriers at the metal-wire contacts. Thus, we conclude 
that the reason for deviation from the maximum con- 
ductance predicted by Landauer theory is backscattering 
from the contact between the electrode and the molecule, 
emphasizing the importance of ab initio studies of inter- 
face chemistry and electron structure. 

In the jellium electrode model, oscillatory behavior can 
be explainedi by the different occupation of molecular 



orbitals. In that model the highest occupied molecu- 
lar orbital (HOMO) in the even-n free chain is half oc- 
cupied, which leads to misalignment of the molecular 
HOMO with the Fermi level of the electrodes, whereas 
the HOMO in odd-n free chains is fully occupied, which 
leads to partial discharge of the molecular HOMO and 
alignment with electrodes Fermi level. Thus, even(odd)- 
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FIG. 5: MPSH states of (a) tetraene and (b) pentaene systems 
around the Fermi level; the corresponding radical orbital is 
referenced according to the level structure on Fig. 2] 
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n chains have low(high) density of states at the Fermi 
level. Comparing with our results with thiolate-capped 
cumulenes chemisorbed onto gold electrodes, we find that 
the molecular level alignment has changed because of 
strong hybridization between molecular and electrode 
surface states coupled by thiolate bonds. The position 
of the LUMO level lies below the Fermi level at zero bias 
for both even and odd numbered cumulene wires. Charge 
transfer from the electrodes to the molecule can further 
partially fill the LUMO. Therefore, the conductance of 
even-n wires is still quite high at zero bias. Comparing 
transmission spectra for odd-n and even-n carbon num- 
ber wires, we see that the main peaks near the Fermi level 
are rather similar, but the peak magnitudes and widths 
are different. 

To elucidate the origin of this behavior, we investi- 
gate how the individual states of cumulene wires con- 
tribute to the conductance. MPSH states on the cumu- 
lenes are therefore analyzed to relate transmission spec- 
tra to molecular orbitals. As indicated by the MPSH 
eigenvalues in Fig. |4l the two main transmission peaks 
(marked 'A' and 'B') of tetraene are mainly contributed 
by states 21 and 22, while the peaks of pentaene take 
their origin from states 22, 23, 24 and 25. The corre- 
sponding MPSH states, which are mainly responsible for 
the transmission around the Fermi level, are shown in 
Fig. [5] It is found that common features of these orbitals 
are their significant delocalization along the backbone of 
the wire. Although the shape of peak A and peak B 
of even-n and odd-n wires are rather similar, the trans- 
mission channels inside the molecular wires are different. 
We find that ir^ states of the LUMO level located at peak 
A and peak B of tetraene have no pronounced appear- 
ance (Fig. [5ji, state 21 and 22), while those of pentaene 
show both TTx and tt^ states (Fig. [SId, states 22, 23, 24 
and 25). Thus the occupation of the tetraene LUMO-l-1 
(state 22) upon adsorption of the molecule between metal 



electrodes becomes unfavorable, due to the close symme- 
try of this state and lower-lying state 21, which keeps the 
state far above Fermi level. The iTy state in tetraene ap- 
pears at MPSH state 20, which is far away from the Fermi 
level. This makes the number of transmission eigenchan- 
nels around the Fermi level of odd-n wires higher than in 
even-n wires. This physical picture of electrode- induced 
interaction and hybridization of states upon chemisorp- 
tion clarifies the origin of the oscillatory conductance in 
cumulenes wires. 

Due to its symmetric structure, the /- V curve and 
conductance of cumulenes are symmetric with bias po- 
larity. As the applied voltage is increased, the conduc- 
tance of cumulenes slowly decreases. This result can be 
explained by plotting transmission spectra as a function 
of bias voltage, as shown in Fig. [B] As positive bias is 
applied, the shapes of transmission spectra remain simi- 
lar, whereas amplitudes of peak A decrease simultaneous 
with the decrease gap between peak A and peak B. The 
change of its position and magnitude relative to the elec- 
trochemical potentials of each electrode results in a grad- 
ual drop of conductance. The obtained curves demon- 
strate metallic IVCs with linear increase in bias range 
from -0.9 to 0.9 V of odd-number wire, and from -0.7 to 
0.7 V of even-number wire. Above those values of bias 
voltage, /- V curves start to deviate from linearity. In 
the linear regime, the slopes of the /- V curves give the 
conductance (G) of the cumulene systems at zero bias, 
which is 1.29 for n = 4, 1.68 for n = 5, 1.42 for n = 6, 
1.59 for n = 7, 1.42 for n = 8, 1.69 for n = 9 in units of 
quantum conductance (Go). Since our electrodes are very 
different from jellium electrodes, the obtained magnitude 
of conductance is slightly lower than obtained in Ref. H. 
Additionally, the conductance change upon stretching of 
the wire was found to be of the order of 0.1 (Gq) for 1 
A change in the inter-electrode distance. 
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FIG. 6: Transmission spectra of pentaene system as a function 
of bias voltage. Energies obtained from the average electro- 
chemical potential of the electrodes. 



IV. SUMMARY 

We have presented calculated IVCs of the cumu- 
lene wire junctions demonstrating ballistic-like transport. 
The high metallic-like conductance varies in an oscilla- 
tory manner with odd-even n: odd-n «1.7 (Go), even- 
n sal. 4 (Go). This study indicates that cumulenes are 
worthy of further theoretical and experimental research, 
with potential to become useful one-dimensional molec- 
ular metallic wires for nano-scale electronic devices. We 
emphasize the decisive role of the anchoring group for 
adsorbed molecular structure and properties. 
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